(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     31737,        990]
NotebookOptionsPosition[     30900,        959]
NotebookOutlinePosition[     31274,        975]
CellTagsIndexPosition[     31231,        972]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Ke", "=", 
   RowBox[{"{", 
    RowBox[{"{", "1", "}"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tl", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], ",", "0", ",", "1", ",", "0"}], "}"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Kl", "=", 
   RowBox[{
    RowBox[{"Tl", "\[Transpose]"}], ".", "Ke", ".", "Tl"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MatrixForm", "[", "Kl", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "x"], ",", 
       SubscriptBox["a", "y"]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "y"]}], ",", 
       SubscriptBox["a", "x"]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"o", "=", 
   RowBox[{"ConstantArray", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{"2", ",", "2"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"ArrayFlatten", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"t", ",", "o"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"o", ",", "t"}], "}"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"K", "=", 
   RowBox[{
    RowBox[{"T", "\[Transpose]"}], ".", "Kl", ".", "T"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MatrixForm", "[", "K", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Nl", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"1", "-", "\[Xi]"}], ",", "0", ",", "\[Xi]", ",", "0"}], "}"}],
      ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"1", "-", "\[Xi]"}], ",", "0", ",", "\[Xi]"}], "}"}]}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Ml", "=", 
   RowBox[{
    RowBox[{"Nl", "\[Transpose]"}], ".", "Nl"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MatrixForm", "[", "Ml", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Ml", "=", 
   RowBox[{"Integrate", "[", 
    RowBox[{"Ml", ",", 
     RowBox[{"{", 
      RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"Ml", "*", "6"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"M", "=", 
   RowBox[{
    RowBox[{"T", "\[Transpose]"}], ".", "Ml", ".", "T"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"Simplify", "[", 
   RowBox[{"M", "*", "6"}], "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.629187086531562*^9, 3.6291871728933125`*^9}, {
   3.6291872229694*^9, 3.6291872309254136`*^9}, {3.629187411979329*^9, 
   3.6291874440217853`*^9}, {3.629187482834653*^9, 3.62918752136672*^9}, {
   3.6291876426101313`*^9, 3.6291876643565693`*^9}, {3.6291877975340014`*^9, 
   3.629187826284852*^9}, {3.6291882982819214`*^9, 3.6291883161929455`*^9}, {
   3.629188347877758*^9, 3.629188401878847*^9}, {3.629188731436696*^9, 
   3.6291887414032664`*^9}, {3.629188784088708*^9, 3.6291888785451107`*^9}, {
   3.629189555729843*^9, 3.629189597529234*^9}, {3.6291896318831987`*^9, 
   3.6291896624509473`*^9}, {3.6291899619960804`*^9, 3.629189993231867*^9}, {
   3.6291900318840775`*^9, 3.629190034275214*^9}, {3.629190240889032*^9, 
   3.629190260022126*^9}, {3.6291903299921284`*^9, 3.629190342944869*^9}, {
   3.629190754775425*^9, 3.6291907665520983`*^9}, {3.629191038238638*^9, 
   3.6291910863123875`*^9}, {3.6291912195180063`*^9, 3.629191223727247*^9}, {
   3.6291912611753893`*^9, 3.629191276082242*^9}, 3.629193417940749*^9, {
   3.629258908501474*^9, 3.6292589088914747`*^9}, {3.62925908390818*^9, 
   3.6292591118322287`*^9}, {3.62926130532382*^9, 3.629261359384912*^9}, {
   3.6292614449178042`*^9, 3.629261471615331*^9}, {3.629262443423916*^9, 
   3.6292624556696157`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", 
      RowBox[{"-", "1"}], "0"},
     {"0", "0", "0", "0"},
     {
      RowBox[{"-", "1"}], "0", "1", "0"},
     {"0", "0", "0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.629188794842323*^9, 3.6291888177866354`*^9}, 
   3.629188880822241*^9, 3.6291896646600733`*^9, 3.629190264241368*^9, 
   3.62919108775247*^9, 3.62919122551635*^9, 3.629191277540325*^9, 
   3.6291934262932267`*^9, {3.629259070570157*^9, 3.62925908975819*^9}, 
   3.629259130739462*^9, 3.629261327478087*^9, 3.629261360447973*^9, 
   3.6292614800408134`*^9, 3.62926244595306*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      SubsuperscriptBox["a", "x", "2"], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "x", "2"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}]},
     {
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      SubsuperscriptBox["a", "y", "2"], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "y", "2"]}]},
     {
      RowBox[{"-", 
       SubsuperscriptBox["a", "x", "2"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      SubsuperscriptBox["a", "x", "2"], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}]},
     {
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "y", "2"]}], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      SubsuperscriptBox["a", "y", "2"]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.629188794842323*^9, 3.6291888177866354`*^9}, 
   3.629188880822241*^9, 3.6291896646600733`*^9, 3.629190264241368*^9, 
   3.62919108775247*^9, 3.62919122551635*^9, 3.629191277540325*^9, 
   3.6291934262932267`*^9, {3.629259070570157*^9, 3.62925908975819*^9}, 
   3.629259130739462*^9, 3.629261327478087*^9, 3.629261360447973*^9, 
   3.6292614800408134`*^9, 3.6292624459590607`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"2", "0", "1", "0"},
     {"0", "2", "0", "1"},
     {"1", "0", "2", "0"},
     {"0", "1", "0", "2"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.629188794842323*^9, 3.6291888177866354`*^9}, 
   3.629188880822241*^9, 3.6291896646600733`*^9, 3.629190264241368*^9, 
   3.62919108775247*^9, 3.62919122551635*^9, 3.629191277540325*^9, 
   3.6291934262932267`*^9, {3.629259070570157*^9, 3.62925908975819*^9}, 
   3.629259130739462*^9, 3.629261327478087*^9, 3.629261360447973*^9, 
   3.6292614800408134`*^9, 3.6292624492322474`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"2", "0", "1", "0"},
     {"0", "2", "0", "1"},
     {"1", "0", "2", "0"},
     {"0", "1", "0", "2"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.629188794842323*^9, 3.6291888177866354`*^9}, 
   3.629188880822241*^9, 3.6291896646600733`*^9, 3.629190264241368*^9, 
   3.62919108775247*^9, 3.62919122551635*^9, 3.629191277540325*^9, 
   3.6291934262932267`*^9, {3.629259070570157*^9, 3.62925908975819*^9}, 
   3.629259130739462*^9, 3.629261327478087*^9, 3.629261360447973*^9, 
   3.6292614800408134`*^9, 3.629262449236248*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Ke", "=", 
   RowBox[{"{", 
    RowBox[{"{", "1", "}"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tl", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], ",", "0", ",", "0", ",", "1", ",", "0", ",", "0"}], 
     "}"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Kl", "=", 
   RowBox[{
    RowBox[{"Tl", "\[Transpose]"}], ".", "Ke", ".", "Tl"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MatrixForm", "[", "Kl", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "x"], ",", 
       SubscriptBox["a", "y"], ",", 
       SubscriptBox["a", "z"]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "x"], ",", 
       SubscriptBox["b", "y"], ",", 
       SubscriptBox["b", "z"]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "x"], ",", 
       SubscriptBox["c", "y"], ",", 
       SubscriptBox["c", "z"]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"o", "=", 
   RowBox[{"ConstantArray", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{"3", ",", "3"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"ArrayFlatten", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"t", ",", "o"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"o", ",", "t"}], "}"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"K", "=", 
   RowBox[{
    RowBox[{"T", "\[Transpose]"}], ".", "Kl", ".", "T"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MatrixForm", "[", "K", "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Nl", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"1", "-", "\[Xi]"}], ",", "0", ",", "0", ",", "\[Xi]", ",", 
       "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"1", "-", "\[Xi]"}], ",", "0", ",", "0", ",", "\[Xi]", ",", 
       "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       RowBox[{"1", "-", "\[Xi]"}], ",", "0", ",", "0", ",", "\[Xi]"}], 
      "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Ml", "=", 
   RowBox[{
    RowBox[{"Nl", "\[Transpose]"}], ".", "Nl"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MatrixForm", "[", "Ml", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Ml", "=", 
   RowBox[{"Integrate", "[", 
    RowBox[{"Ml", ",", 
     RowBox[{"{", 
      RowBox[{"\[Xi]", ",", "0", ",", "1"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"Ml", "*", "6"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"M", "=", 
   RowBox[{
    RowBox[{"T", "\[Transpose]"}], ".", "Ml", ".", "T"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"Simplify", "[", 
   RowBox[{"M", "*", "6"}], "]"}], "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.6292589231342993`*^9, 3.6292589260203047`*^9}, {
  3.62925929628695*^9, 3.629259297534952*^9}, {3.6292593301234093`*^9, 
  3.629259377469492*^9}, {3.629266633531246*^9, 3.6292666799569273`*^9}, {
  3.6292667119525833`*^9, 3.6292667208601985`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0", 
      RowBox[{"-", "1"}], "0", "0"},
     {"0", "0", "0", "0", "0", "0"},
     {"0", "0", "0", "0", "0", "0"},
     {
      RowBox[{"-", "1"}], "0", "0", "1", "0", "0"},
     {"0", "0", "0", "0", "0", "0"},
     {"0", "0", "0", "0", "0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6292589313399143`*^9, 3.6292593842087035`*^9, 
  3.6292667348690233`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      SubsuperscriptBox["a", "x", "2"], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "x", "2"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "z"]}]},
     {
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      SubsuperscriptBox["a", "y", "2"], 
      RowBox[{
       SubscriptBox["a", "y"], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "y", "2"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "y"]}], " ", 
       SubscriptBox["a", "z"]}]},
     {
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{
       SubscriptBox["a", "y"], " ", 
       SubscriptBox["a", "z"]}], 
      SubsuperscriptBox["a", "z", "2"], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "y"]}], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "z", "2"]}]},
     {
      RowBox[{"-", 
       SubsuperscriptBox["a", "x", "2"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "z"]}], 
      SubsuperscriptBox["a", "x", "2"], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "z"]}]},
     {
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "y"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "y", "2"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "y"]}], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "y"]}], 
      SubsuperscriptBox["a", "y", "2"], 
      RowBox[{
       SubscriptBox["a", "y"], " ", 
       SubscriptBox["a", "z"]}]},
     {
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "x"]}], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{
       RowBox[{"-", 
        SubscriptBox["a", "y"]}], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{"-", 
       SubsuperscriptBox["a", "z", "2"]}], 
      RowBox[{
       SubscriptBox["a", "x"], " ", 
       SubscriptBox["a", "z"]}], 
      RowBox[{
       SubscriptBox["a", "y"], " ", 
       SubscriptBox["a", "z"]}], 
      SubsuperscriptBox["a", "z", "2"]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6292589313399143`*^9, 3.6292593842087035`*^9, 
  3.629266734915823*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"2", "0", "0", "1", "0", "0"},
     {"0", "2", "0", "0", "1", "0"},
     {"0", "0", "2", "0", "0", "1"},
     {"1", "0", "0", "2", "0", "0"},
     {"0", "1", "0", "0", "2", "0"},
     {"0", "0", "1", "0", "0", "2"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6292589313399143`*^9, 3.6292593842087035`*^9, 
  3.629266735664624*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SubsuperscriptBox["a", "x", "2"], "+", 
         SubsuperscriptBox["b", "x", "2"], "+", 
         SubsuperscriptBox["c", "x", "2"]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "y"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "y"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "y"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}], 
      RowBox[{
       SubsuperscriptBox["a", "x", "2"], "+", 
       SubsuperscriptBox["b", "x", "2"], "+", 
       SubsuperscriptBox["c", "x", "2"]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "y"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "y"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "y"]}]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "z"]}]}]},
     {
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "y"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "y"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "y"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SubsuperscriptBox["a", "y", "2"], "+", 
         SubsuperscriptBox["b", "y", "2"], "+", 
         SubsuperscriptBox["c", "y", "2"]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "y"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "y"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "y"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "y"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "y"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "y"]}]}], 
      RowBox[{
       SubsuperscriptBox["a", "y", "2"], "+", 
       SubsuperscriptBox["b", "y", "2"], "+", 
       SubsuperscriptBox["c", "y", "2"]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "y"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "y"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "y"], " ", 
        SubscriptBox["c", "z"]}]}]},
     {
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "y"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "y"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "y"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SubsuperscriptBox["a", "z", "2"], "+", 
         SubsuperscriptBox["b", "z", "2"], "+", 
         SubsuperscriptBox["c", "z", "2"]}], ")"}]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "z"]}]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "y"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "y"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "y"], " ", 
        SubscriptBox["c", "z"]}]}], 
      RowBox[{
       SubsuperscriptBox["a", "z", "2"], "+", 
       SubsuperscriptBox["b", "z", "2"], "+", 
       SubsuperscriptBox["c", "z", "2"]}]},
     {
      RowBox[{
       SubsuperscriptBox["a", "x", "2"], "+", 
       SubsuperscriptBox["b", "x", "2"], "+", 
       SubsuperscriptBox["c", "x", "2"]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "y"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "y"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "y"]}]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "z"]}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SubsuperscriptBox["a", "x", "2"], "+", 
         SubsuperscriptBox["b", "x", "2"], "+", 
         SubsuperscriptBox["c", "x", "2"]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "y"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "y"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "y"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}]},
     {
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "y"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "y"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "y"]}]}], 
      RowBox[{
       SubsuperscriptBox["a", "y", "2"], "+", 
       SubsuperscriptBox["b", "y", "2"], "+", 
       SubsuperscriptBox["c", "y", "2"]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "y"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "y"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "y"], " ", 
        SubscriptBox["c", "z"]}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "y"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "y"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "y"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SubsuperscriptBox["a", "y", "2"], "+", 
         SubsuperscriptBox["b", "y", "2"], "+", 
         SubsuperscriptBox["c", "y", "2"]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "y"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "y"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "y"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}]},
     {
      RowBox[{
       RowBox[{
        SubscriptBox["a", "x"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "x"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "x"], " ", 
        SubscriptBox["c", "z"]}]}], 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "y"], " ", 
        SubscriptBox["a", "z"]}], "+", 
       RowBox[{
        SubscriptBox["b", "y"], " ", 
        SubscriptBox["b", "z"]}], "+", 
       RowBox[{
        SubscriptBox["c", "y"], " ", 
        SubscriptBox["c", "z"]}]}], 
      RowBox[{
       SubsuperscriptBox["a", "z", "2"], "+", 
       SubsuperscriptBox["b", "z", "2"], "+", 
       SubsuperscriptBox["c", "z", "2"]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "x"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "x"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "x"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "y"], " ", 
          SubscriptBox["a", "z"]}], "+", 
         RowBox[{
          SubscriptBox["b", "y"], " ", 
          SubscriptBox["b", "z"]}], "+", 
         RowBox[{
          SubscriptBox["c", "y"], " ", 
          SubscriptBox["c", "z"]}]}], ")"}]}], 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SubsuperscriptBox["a", "z", "2"], "+", 
         SubsuperscriptBox["b", "z", "2"], "+", 
         SubsuperscriptBox["c", "z", "2"]}], ")"}]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6292589313399143`*^9, 3.6292593842087035`*^9, 
  3.6292667357114244`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"IdentityMatrix", "[", "6", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.629268150477891*^9, 3.6292681947663684`*^9}, {
  3.6292682354824395`*^9, 3.6292682499904647`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0", "0", "0", "0"},
     {"0", "1", "0", "0", "0", "0"},
     {"0", "0", "1", "0", "0", "0"},
     {"0", "0", "0", "1", "0", "0"},
     {"0", "0", "0", "0", "1", "0"},
     {"0", "0", "0", "0", "0", "1"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.629268251956068*^9}]
}, Open  ]]
},
WindowSize->{1350, 709},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
Magnification:>2. Inherited,
FrontEndVersion->"10.0 for Microsoft Windows (64-bit) (September 9, 2014)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 4019, 111, 707, "Input"],
Cell[4602, 135, 1103, 27, 146, "Output"],
Cell[5708, 164, 2186, 65, 174, "Output"],
Cell[7897, 231, 1062, 25, 146, "Output"],
Cell[8962, 258, 1060, 25, 146, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10059, 288, 3429, 107, 745, "Input"],
Cell[13491, 397, 912, 25, 233, "Output"],
Cell[14406, 424, 3668, 125, 273, "Output"],
Cell[18077, 551, 866, 23, 233, "Output"],
Cell[18946, 576, 10858, 347, 273, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29841, 928, 225, 4, 60, "Input"],
Cell[30069, 934, 815, 22, 290, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
